Constraining invisible neutrino decays with the cosmic microwave background 
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Precision measurements of the acoustic peaks of the cosmic microwave background indicate that 
neutrinos must be freely streaming at the photon decoupling epoch when T « 0.3 eV. This re- 
quirement implies restrictive limits on "secret neutrino interactions," notably on neutrino Yukawa 
couplings with hypothetical low-mass (pseudo)scalars (j>. For diagonal couplings in the neutrino 
mass basis we find g < f x 10 -7 , comparable to limits from supernova 1987A. For the off-diagonal 
couplings and assuming hierarchical neutrino masses we find g < f x 10 -11 (0.05 eV/m) 2 where m 
is the heavier mass of a given neutrino pair connected by g. This stringent limit excludes that the 
flavor content of high-energy neutrinos from cosmic-ray sources is modified by v — ► v' + decays 
on their way to Earth. 



PACS numbers: 14.60.Lm, 98.80.-k, 95.35.+d 



I. INTRODUCTION 

The observed acoustic peaks in the temperature distri- 
bution of the cosmic microwave background (CMB) by 
WMAP pi Q and other experiments such as CBI 
DASI @, ACBAR and BOOMERANG H00 h &ve 
provided a plethora of detailed information about our 
universe. In particular, several authors have indepen- 
dently realized that these observations imply that neu- 
trinos must be freely streamin g ar ound the photon decou- 
pling epoch at T » 0.3 eV (the fact that neutrino 
free-streaming affects the CMB acoustic peaks in a very 
characteristic way was first discussed in [llllj). While 
ordinary weak interactions freeze out at T w 1 MeV, neu- 
trinos could have "secret interactions" that are still in 
equilibrium at late times or that actually recouple at late 
times. In particular, this applies to neutrino interactions 
with new massless or low-mass scalars or pseudoscalars. 
Typically, these particles would be the Nambu-Goldstone 
bosons of a new symmetry that is broken at some low 
energy scale. The Majoron model and its variants is a 
widely discussed case in point 0, 0, 0, 0] . 

The late free-streaming requirement allows one to test 
neutrino interactions at eV energies, far below what 
is possible in the laboratory. It was previously recog- 
nized that signatures for neutrino mass generation in 
Majoron type models may show up in future CMB ob- 
servations |l2j and conversely, that existing observations 
exclude the "neutrinoless universe" that had been in- 
voked to escape the cosmological neutrino mass limit |l8j . 

The purpose of our paper is to show that the free- 
streaming requirement translates into very stringent lim- 
its on the neutrino-(pseudo)scalar Yukawa couplings. 
Our limits suggest that interactions of this sort play 



no significant role for supernova physics 0, |2£l l2l| . 

Perhaps more interestingly, our limits exclude scenar- 
ios where the flavor content of high-energy cosmic-ray 
neutrinos is modified by decays v — > v' + (f> in addi- 
tion to the standard modification caused by flavor os- 
cillations ihei. 

In Sec. II we first consider binary interactions among 
neutrinos and bosons that lead to new limits on the di- 
agonal and off-diagonal Yukawa couplings. In Sec. Ill 
we consider decay and coalescence processes of the type 
1^2 + 3 that provide new and very restrictive limits on 
the off-diagonal interactions. We summarize our findings 
in Sec. IV. 



II. BINARY INTERACTIONS 

We study ordinary neutrinos interacting with low-mass 
pseudoscalars <j) via the coupling 



L 



jk 



(1) 
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If the pseudoscalars are Nambu-Goldstone bosons of a 
new symmetry, as one should expect, then a derivative 
coupling is more appropriate, but for our most interesting 
process, neutrino decay v — ► v' + tf>, the pseudoscalar and 
derivative couplings are equivalent. For binary processes 
it turns out to be conservative to use the pseudoscalar 
coupling, a point discussed in more detail at the end of 
this section. We do not explicitly study scalar interac- 
tions, but the results would be quantitatively similar. 

The interaction Eq. allows for binary processes of 
the type v + D <-> <fi + <fi, v + 4> <-> u + <f>, and finally 
v + v — > v + v with a 0-exchange. For this latter process 
the pseudoscalar and derivative couplings are equivalent 
because each fermion line has only one Nambu-Goldstone 
boson attached to it. Apart from numerical factors the 
scattering rate in a thermal environment of relativistic 
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neutrinos is 

rW. 3 +4~.9 4 T, (2) 

where g is the largest entry of the Yukawa coupling 
matrix. To avoid acoustic oscillations of the neutrino- 
Majoron fluid we require that the cosmic expansion 
rate H exceeds I\ + 2^3+4 at the time of photon decou- 
pling. Of course, this criterion is somewhat schematic, 
but to be numerically specific we will use the follow- 
ing conditions. The photon temperature at decou- 
pling is T 7i dcc = 0.256 eV. The corresponding neu- 
trino temperature relevant for our estimate is T„ ; dcc = 
(4/ll) 1 / 3 T 7jd ec = 0.18 cV. At this epoch the uni- 
verse is matter dominated so that the expansion rate is 
H dcc = 100 km s- 1 Mpc" 1 (Q M h 2 ) 1/2 (zd C c + 1) 3/2 - With 
Om/i 2 = 0.134 and z dcc = 1088 for the cosmic matter 
density and the redshift at decoupling, respectively Q, 
we find 

iJdcc = 4.27 x 10~ 14 s" 1 = 2.81 x 10~ 29 eV . (3) 

The free-streaming requirement T 1+2 ^3+4 < -Hdec thus 
translates into 

5<l.lxl0" 7 . (4) 

Because of the g 4 dependence of the interaction rate, cor- 
rections to this limit from exact numerical factors are 
minor, and the limit correspondingly robust. 

Our limit Eq. @ on the largest of the Yukawa cou- 
plings is nominally more restrictive than the limit ob- 
tained from the energy-loss argument of supernova (SN) 
1987A |l9ll2Ct l2l|. However, since our limit is based on 
a dimensional analysis with uncertain numerical factors, 
one can only claim that the two limits are comparable. 

If the bosons are Nambu-Goldstone bosons of a new 
broken symmetry, there is a relation g = m/f between 
the diagonal Yukawa couplings, the neutrino masses, 
and the symmetry breaking scale /. For processes 
with two Nambu-Goldstone boson lines attached to one 
fermion line, the coupling is derivative rather than pseu- 
doscalar. The relevant cross section will be proportional 
to m 2 E 2 / f 4 with E a typical energy of the process 
rather than proportional to g 4 — m 4 / f 4 that we used 
in our estimate. If neutrinos are relativistic at T^dec, 
evidently the correct derivative structure for the inter- 
action would lead to more restrictive limits. We note 
that current cosmological limits on the neutrino masses 
are yi m < 1—1.5 eV, and even J2m < 0.4 eV has been 
claimed |22j . Typical neutrino energies at photon decou- 
pling are E « 3T„ ircc « 1 eV so that even in the case of 
degenerate neutrino masses it is justified to treat neutri- 
nos as relativistic. Note that the mass bounds have been 
obtained assuming non-interacting neutrinos, but since 
strongly interacting neutrinos are excluded by CMB the 
bound is self-consistent. 

Of course, if the boson masses were much larger than 
the eV-scale, our considerations would change and these 
restrictive limits could be avoided. 



III. DECAY AND COALESCENCE 

Of greater interest is the decay process v — ► v' + <f> and 
the coalescence process v 1 + <fi — * v as these are only of 
first order and thus the rate is proportional to g 2 rather 
than g 4 . Here it is exact to use the non-derivative form of 
the interaction. The sum of the decay rates for v — » v'+(j) 
and v — > v 1 + cp in the rest frame of the parent neutrino 
with mass m m' is [2^, |26| 

. 2 

Tdecay = TTl ■ (5) 

In the frame of the thermal medium, a typical neutrino 
energy is E = 3T so that the rate is reduced by a typical 
Lorentz factor m/3T. 

The decay and coalescence processes are kinematically 
constrained, for relativistic particles, to couple nearly 
collinear modes of the interacting particles. Therefore, 
even if the decay is isotropic in the rest frame of the 
parent particle, the decay products will have directions 
within an approximate angle corresponding to the par- 
ent's Lorentz factor m/E. Therefore, to randomize the 
direction of the original neutrino requires a random walk 
of small angular steps. Therefore, we must include 
a factor {m/E) 2 in the medium- frame interaction rate 
to obtain the relevant "transport rate" rather than a 
naive interaction rate. The same argument was raised 
in Ref. 01 ■ We also note that the usual transport cross 
section is the ordinary cross section times (1 — cos 8) with 
8 the scattering angle. Therefore, for small angles the 
transport cross section is the ordinary cross section times 
\8 2 , resulting in the same approximate (m/E) 2 correc- 
tion factor. 

We conclude that we should compare the "transport 
rate" 

with the expansion rate at photon decoupling Eq. J2J. 
Using E = ST^dec we thus find 

S < 12(3»)" J (tf dot r,J d „ c ) 1/2 m- J 

= 0. 6 lxlO-"(^)\ (7) 

where as a mass scale we have used the largest neutrino 
mass of about 50 meV implied by oscillation data in a 
hierarchical mass scenario. This is by far the most re- 
strictive limit on the off-diagonal neutrino-Majoron cou- 
plings. 

Translating this constraint into a limit on the rest- 
frame lifetime yields 

^ 2 > <10 '° * (m^v) 3 ■ < 8 » 

Therefore, neutrinos could still be rather short-lived on 
cosmological time scales. 
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If the coupling g is between two nearly degenerate mass 
eigenstates, the decay rate acquires an additional factor 
of approximately (<5m 2 ) 3 /m 6 [2(| so that the limit on g 
is relaxed by a factor m 3 /(<5m 2 ) 3 / 2 . 

IV. CONCLUSIONS 

We have shown that the neutrino free-streaming re- 
quirement at the photon decoupling epoch implies new 
limits on the neutrino Yukawa couplings with hypotheti- 
cal low-mass bosons. For the diagonal couplings we find 
g < 10~ 7 , comparable to the energy-loss argument of 
SN 1987A [H Hi E3- The SN limits suffer from the 
complication that one needs to distinguish between the 
free-streaming regime where Majoron emission simply 
provides a new channel of energy loss and the trapping 
regime where neutrino-Majoron interactions are so strong 
that Majorons are trapped. Unless the Majorons or other 
new bosons have masses which are large compared to 
the eV scale, our limits imply that Majorons and simi- 
lar particles interact too weakly to be trapped in a SN 
core. Given the uncertainty of our limit, Majoron emis- 
sion could perhaps still provide a non-negligible channel 

I 



> 1.1 x 10- ( *™*) 
\ m J 



Therefore, g would need to be about four orders of mag- 
nitude larger than our new limit. Even if we consider 
decays between the second-lightest and lightest neutri- 
nos with m ss 10 meV and/or larger distances or smaller 
energies, this conclusion is not changed. For degenerate 
neutrino masses, the decay rate of cosmic-ray neutrinos 
and the early-universe reaction rate get penalized with 
the same factor so that, again, our conclusion remains 
unchanged. 

Therefore, it appears that possible future deviations 
from the expected 1:1:1 flavor content of high-energy 
cosmic-ray neutrinos would have to be ascribed to other 
effects than Majoron decays. By the same token, invisible 
neutrino decays of the Majoron type can not affect solar 
or atmospheric neutrino observations. 

In essence the reason why our bound is so much 
stronger is that neutrinos are almost non-relativistic at 
decoupling, having energies in the sub-eV range. The 
Lorentz factor is therefore enormously smaller than the 
one considered for TeV neutrinos. On the other hand the 
effective "baseline" for decays is roughly H^ ec ~ 0.2 Mpc, 
a number which is only 3 orders of magnitude smaller 
than the 100 Mpc considered typical for high energy neu- 
trinos observable in neutrino telescopes. 

The diffuse cosmic neutrino background from all core- 



of energy loss, but a dominant dynamical role appears to 
be excluded. 

We obtain a much more restrictive limit on the off- 
diagonal couplings of g < 10~ n (50 meV/m) 2 with m 
the heavier mass of a given pair of neutrinos with non- 
degenerate masses. This is by far the most restrictive 
limit on such interactions. 

It is of interest in the context of scenarios where high- 
energy neutrinos from cosmic-ray sources may decay on 
their way to Earth. These neutrinos are produced by 
the decay of charged pions that in turn are produced as 
secondary products by the high-energy cosmic-ray pro- 
tons that must be produced somewhere in the universe. 
Standard flavor oscillations imply an expected flavor ra- 
tio at Earth of u e : : v r « 1 : 1 : 1. Therefore, 
if future large-scale neutrino telescopes were to observe 
significant deviations from this expected flavor composi- 
tion one would be tempted to conclude that these mod- 
ifications are caused by v — > v' + <fi decays [23|, ■ As- 
suming a source distance of D = 100 Mpc and a neu- 
trino energy E = 10 TeV, a strong decay effect obtains 
if T decay (m/E) > D~ x . With Eq. © this implies the 
requirement 



i 

collapse supernovae in the universe is another potentially 
detectable neutrino flux from a cosmological distance. 
This flux can also be affected by invisible neutrino de- 
cays [27], 13 ■ Typical energies are in the range of tens of 
MeV, i.e. six orders of magnitude smaller relative to our 
discussion of high-energy neutrinos. The typical source 
distance is at least a factor of ten larger. Overall the 
required coupling strength for significant decay effects to 
occur is 3-4 orders of magnitude smaller than Eq. 10. 
Therefore, given the crude nature of our limit we can not 
claim with complete confidence that the diffuse super- 
nova neutrinos are not affected by Majoron-type decays. 

Either way, the universe as a neutrino laboratory once 
more provides information on these no-longer-so-elusive 
particles which is of direct relevance to other experimen- 
tal directions in neutrino research. 
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